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Abstract 

We report on searches for ^^'^ ^ ~^ Dg^cj). In the context of the Standard Model 

the branching fractions for these decays are expected to be highly suppressed, since they proceed 
through annhilation of the b and u quarks in the B~ meson. Our results are based on 234 million 
T(42) BE decays collected with the BABAR detector at SLAG. We find no evidence for these 
decays, and we set 90% confidence level upper limits on the branching fractions 

B{B- -> Dg^) < 1.8 X 10-6 

B{B- D*f^) < 1.1 X 10^^ 
These results are consistent with Standard Model expectations. 
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1 Introduction 



In the Standard Model (SM), the decay B~ — > Dg (j) occurs through annhilation of the two 
quarks in the 5-meson into a virtual W , see Figure 1. No annhilation- type B decays have ever 
been observed to date. The current 90% C.L. upper limits on B~ and B~ (j) are 

3.2 X 10-"^ and 4 x IQ-^, respectively [1]. 



B- 




Figure 1: Feynman diagram for B Dg (j). 



In the SM, annhilation diagrams are highly suppressed. Calculations of the B^ branch- 
ing fraction give predictions of 3 x 10^'^ using a perturbative QCD approach [2], or 7 x 10^^ using 
QCD-improved factorization [3]. 

Since the current experimental limits are about three orders of magnitude higher than the SM 

— (*)— 

expectations, searches for B — > (j) could be sensitive to new physics contributions. For 

example, in Reference [3] the branching fraction of B~ D^cp is estimated to be 8 x 10^^ in a two 
Higgs doublet model and 3 x 10"'^ in the minimal supersymmetric model with i?-parity violation. 

2 The BaBAR detector and dataset 

Our results are based on 234 x 10^ T(4S') BB decays, corresponding to an integrated luminosity 
of 212 fb-^ cohected between 1999 and 2004 with the BABAR detector [4] at the PEP-II B Factory 
at SLAC [5]. A 12 fb~^ off-resonance data sample, with a center of mass (CM) energy 40 MeV 
below the T(4<S') resonance peak, is used to study continuum events, e~^e~ qq (q = u,d,s, or 
c). The number of B-mesons in our data sample is two orders of magnitudes larger than in the 
previously published search for B~ D^g ^ (p [1]. 

3 Analysis method 

— (*)— *— — — — — 

We search for the decay B Dg (f) in the following modes: Dg ^5 7) ~^ ^'^ ' , 
and K*^K~, (p — > K^K~, Ks vr^vr", and K*^ K^tt^ (charged conjugate decay modes are 
implied throughout this article). We denote the (p from B^ (p decay as the "bachelor (p" , 

in order to distinguish it from the (p in the D'^ — > (pTT~ decay. 
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All kaon candidate tracks in the reconstructed decay chains must satisfy a set of loose kaon 
identification criteria based on the response of the internally-reflecting ring-imaging Cherenkov 
radiation detector and the ionization measurements in the Drift Chamber and the Silicon Vertex 
Tracker. The kaon selection efficiency is a function of momentum and polar angle, and is typically 
95%. These requirements provide a rejection factor of order 10 against pion backgrounds. No 
particle identification requirements are imposed on pion candidate tracks. 

We select (f), Kg, and candidates from pairs of oppositely-charged tracks with invariant 
masses consistent with the parent particle decay hypothesis and consistent with originating from a 
common vertex. The invariant mass requirements are ±10 MeV {j^ 2.4r) for the ^, ±9 MeV (~ Sa) 
for the Ks, and ±75 MeV (~ 1.5r) for the K*^ . We then form Dg candidates in the three modes 
listed above by combining ^, Kg, or candidates with an additional track. The invariant mass 
of the candidate must be within 15 MeV (~ 3o") of the known mass. In the (f)TT~ and 

Dg K*^K~ modes, all three charged tracks are required to originate from a common vertex. In 
the KsK~ mode, the Ks and Dg vertices are required to be separated by at least 3 mm. 

This last requirement is very effective in rejecting combinatorial background and is 94% efficient 
for signal. We select D*^ candidates from and photon candidates. The photon candidates are 
constructed from calorimeter clusters with lateral profiles consistent with photon showers and with 
energy above 60 MeV in the laboratory frame. We require that the mass difference AM between 
the D*^ and Dg candidates be between 130 and 156 MeV. The AM resolution is about 5 MeV. 

At each stage in the reconstruction chain, the measurement of the momentum vector of an 
intermediate particle is improved by refitting the momenta of the decay products with kinematical 
constraints. These constraints are based on the known mass of the intermediate particle and on 

the fact that the decay products must originate from a common point in space. 

_ (*)— _ 

Finally, we select B~ candidates by combining Dg ' and bachelor ^ candidates, k candi- 
date is characterized kinematically by the energy-substituted mass rriES = \j {\s + Po ■ pb)'^/Eq — p% 
and energy difference AE = — ^^/s, where E and p are energy and momentum, the asterisk 
denotes the CM frame, the subscripts and B refer to the initial T{4S) and B candidate, respec- 
tively, and s is the square of the CM energy. In the CM frame, ttt-es reduces to tbes = 'sj I ~ P*b ■ 
For signal events we expect mES ~ Mb, the known B" mass, and AE ~ 0. The resolutions on 
ruES and AE are approximately 2.6 MeV and 10 MeV, respectively. 

If there is more than one B~ candidate in an event, we retain the best candidate based on a 
algorithm that uses the measured values, known values, and resolutions for the Dg mass, the 
bachelor ^ mass, and, where applicable, AM. 

Studies of simulated events and off-resonance data indicate that most of the backgrounds to 

— (*)— 

the B — >■ Dg ' signal are from continuum events. To reduce these backgrounds we make two 
additional requirements. First, we require | cos^tI < 0.9, where is the angle between the thrust 
axes of the B~ candidate and the rest of the tracks and neutral clusters in the event, calculated 
in the CM frame. The distribution of | cos Qt\ is essentially uniform for signal events and strongly 
peaked near one for continuum events. Second, for each event we define a relative likelihood 
for signal and background based on a number of kinematical quantities. The relative likelihood is 
defined as the ratio of the likelihoods for signal and background. The signal (background) likelihood 
is defined as the product of the probability density functions, PDFs, for the various kinematical 
quantities in signal (background) events. 

The kinematical quantities used in the likelihood are reconstructed masses, decay angles, and 
a Fisher discriminant designed to distinguish between continuum and BB events. All PDFs are 
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chosen based on studies of Monte Carlo and off-resonance data. 

The masses used in the hkeUhoods are those of the Dg, the (f) in the Dg 4>7r~ , the 
in Dg K*^K~, and AM in Dg~ Dsj. The signal PDFs for the mass variables are the 
sum of two Gaussian distributions for Ds and AM, a Breit Wigncr distribution for the K*^, and 
a Voigtian distribution for the cp from Ds decay. The background PDFs are constant functions. 
Note that the mass of the bachelor ^ and the mass of the Ks in Dg KsK~ are not used in 
the definition of the likelihoods. This is because studies of background event samples suggest that 
background events contain mostly real bachelor (f) and real Kg mesons. 

The decay angles used in the likelihood are those in the K*^ K~^n~ and in the — > K~^K~ 
decay, both for bachelor ^ and the cf) from Dg 07r~ decay. The signal PDFs for these quantities 
are set by angular momentum conservation to be proportional to cos^ 0, where 9 is the decay angle 
for the process. The one exception is the decay angle distribution of the bachelor (j) in B~ D*g~ (j) 
decay, where the polarization of the two vector mesons in the final state is not known. For this 
reason, the decay angle of the bachelor (f) is not used in the definition of the likelihood for the 
B~ — > D*^(f) mode. The background PDFs for these variables are constant functions of cos^. In 
addition, in the likelihood we also use the polar angle of the B~ candidate in the CM frame {Ob)- 
The signal is expected to follow a sin^ 6b distribution, while the background is independent of 
cos Ob- 

The final component of the likelihood is a Fisher discriminant constructed from the quantities 

Lq = X^i Vi and L2 = X^i Pi cos^ Oi calculated in the CM frame. Here, pi is the momentum and Oi 
is the angle with respect to the thrust axis of the B~ candidate of tracks and clusters not used to 
reconstruct the B~ . The signal and background PDFs for this variable are modelled as bifurcated 
Gaussians with different means and standard deviations. Note that this Fisher discriminant is 
highly correlated with the | cos Ot\ variable defined above. It is because of this correlation that the 
jcos^T'l variable is treated separately and not included in the likelihood. 

The combined efficiency of the requirements on likelihood and |cos^r| varies between 71 and 
83%, depending on the mode. These requirements provide a rejection factor of about 7 against 
backgrounds. They were chosen from studies of off-resonance data as well as simulated background 
and signal events. 

After applying the requirements on relative likelihood and | cos Ot\, we also demand that /S.E be 
within 30 MeV (~ 3(t) of its expected mean value for signal events. This mean value is determined 
from simulation, and varies between —3 and MeV, depending on the mode. 

Table 1: Efficiency (ej), branching fractions (BRj), and products of efficiency and branching frac- 

— (*)— 

tions for the modes used in the B Dg (f> search. The uncertainties on the and BRj are 
discussed in the text. Here BRj is the product of branching fractions for the secondary and tertiary 
decays in the i-th decay mode. 



Mode 






ej 


BRi 






BRj 




B- 


Dg^, 


Dg (pTT 


0.192 


11.6 • 10" 


-3 


2.22 


• 10- 


-8 


B- 


Ds<t>, 


Dg K'Ks 


0.177 


8.20 • 10- 


-3 


1.45 


• 10- 


-3 


B- 


Dg^, 


Dg K*^K- 


0.140 


14.5 • 10" 


-3 


2.03 


• 10- 


-3 


B- 


D*g-^, 


Dg <pTT~ 


0.109 


10.9 • 10" 


-3 


1.19 


• 10- 


-3 


B- 


D*g-<P, 


Dg K-Ks 


0.100 


7.70 • 10" 


-3 


0.77 


• 10- 


-3 


B- 


D*g^4>, 


Dg K*^K- 


0.083 


13.6 ■ 10" 


-3 


1.14 


■ 10" 


-3 
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The efficiencies of our selection requirements, shown in Table 1, are determined from simulations. 
In the case of the B~ — > D'^'cj) mode we take the average of the efficiencies calculated assuming 
fully longitudinal or transverse polarization for the two vector meson final state. These efficiencies 
arc found to be the same to within 1%. The quantities BRj in Table 1 are the product of the known 
branching fractions for the secondary decay modes. These are taken from the compilation of the 
Particle Data Group [6] , with the exception of the branching fraction for the Ds (pir mode, for 
which we use the latest most precise measurement B{Ds — ^tt) = (4.8 it 0.6)% [7]. Since the 
branching fractions for the other two Ds modes are measured with respect to the Dg (j)7r mode, 
we have also rescaled their tabulated values from the Particle Data Group accordingly. 

4 Systematic studies 

The systematic uncertainties on the products of efficiency and branching ratio for the secondary 
decays in the decay chain of interest are summarized in Table 2. The largest systematic uncertainty 
is associated with the uncertainty on the Ds 07r branching ratio, which is only known to 12% [7], 
and which is used to normalize all other Ds branching ratios. 

From a purely experimental point of view, the most important uncertainty is due to the uncer- 
tainty in the efficiency of the kaon identification requirements. The efficiency of these requirements 
is calibrated using a sample of kinematically identified D*^ — D^tt"*", D° — K~Tr~^, and is known 
at the level of 2%. Thus, this uncertainty result in a systematic uncertainty of 8% for the efficiency 
of the modes with four charged kaons, i.e. those with Dg K*^K~ and D~^ (j)TT, and 6% for the 
mode with three charged kaons {Dg — > KsK~). A second class of uncertainties is associated with 
the detection efficiency for tracks and clusters in the BABAR detector. Prom studies of a variety of 
control samples, the tracking efficiency is understood at the level of 1.4% or 0.6% for transverse 
momenta below or above 200 MeV/c. There is also a 1.9% uncertainty associated wih the recon- 
struction of the Ks TT~^TT~ decay which can occur a few centimeters away from the interaction 
point. Given the multiplicity and momentum spectrum of tracks in the decay modes of interest, the 
uncertainty on the efficiency to reconstruct the tracks in the 5-decay chain is estimated to be 3.7%. 
In the B~ — > Dg~(f) search there is an additional uncertainty of 1.8% due to the uncertainty on the 
efficiency to reconstruct the photon in the Ds^ decay, and also a 1% uncertainty from the 

unknown polarization in the final state. Finally, to ascertain the systematic due to the efficiency 
of the other event selection requirements, we compute the following efficiency variations: shifting 
the AE by 3MeV (0.3%); shifting the mean of the Ds and 4> masses and AM by 1 MeV (0.2%, 
0.1%, 0.2%, respectively); increasing the width of the Ds and (p masses and AM by 1 MeV (1.5%, 
0.4%, 1.5%, respectively); using a Fisher distribution obtained from the data sample of a similar 
analysis, B — Dtt with D — Kir (3%). Thus we assign a 5% systematic on the combined efficiency 
of these selection criteria. 

5 Physics results 

We determine the yield of signal events from an unbinned extended maximum-likelihood fit to 
the w-Es distribution of B~ candidates satisfying all of the requirements listed above. We fit 
simultaneously in two \AE\ regions: In the signal region the distribution is parametrized as a 
Gaussian and the combinatorial background as a threshold function [8]; in a sideband of AE 
{\AE\ < 200 MeV, excluding the signal region) we fit solely for the ( parameter of the threshold 
function. In our fit, the amplitude of the Gaussian is allowed to fluctuate to negative values, but. 
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Table 2: Systematic uncertainties on J2i ' BRj, where the index i runs over the three Dg modes 
used in this analysis, are the experimental efficiencies and BRj are the branching fractions for 
the i-th mode. 





B- Dad) 


B- DT6 


Ds branching fraction 


14% 


14% 


Dg branching fraction 




2.5% 


Other branching fractions 


1.5% 


1.5% 


Charged kaon ID 


7.5% 


7.5% 


Selection requirements 


5% 


5% 


Tracking and Ks efficiency 


3.7% 


3.7% 


Photon efficiency 




1.8% 


Final state polarization 




1% 


Simulation statistics 


0.6% 


0.6% 


Total 


17% 


17% 



for reasons of numerical stability, the sum of the Gaussian and the threshold function is constrained 
to be positive over the full rriES fit range. The mean and the standard deviations of the Gaussian 
are constrained to the values determined from Monte Carlo simulation. The fitting procedure was 
extensively tested with sets of simulated data, and was found to provide an unbiased estimate of 
the signal yield. 

Figure 2 shows the ttt-es distribution of the selected candidates. We see no evidence for B~ 
D^g^"^. The fitted event yields are N = -1.61^;^ and N = SAtH for the B' Dgcj) and 
B~ ^*s~^ modes, respectively, where the quoted uncertainties correspond to changes of ^ in the 
log-likelihood for the fit. The likelihood curves are shown in Figure 3. The requirement that the 
sum of the Gaussian and the threshold function be always positive results in an effective constraint 
N > —1.6 in the B Dgcf) mode. This is the source of the sharp edge at iV = — 1.6 in the 
likelihood distribution of Figure 3(a). 

We use a Bayesian approach with a flat prior to set 90% confldence level upper limits on the 
branching fractions for the B~ — > Dg(j) and B~ modes. In a given mode, the upper limit 

on the number of observed events (Nm) is defined as 

I jr{N)dN = jr{N)dN (1) 

where C{N) is the likelihood as a function of the number of signal events N as determined from 
the mES fit described above. Then the upper limit B on the branching fraction is 

^BB J2i ■ BRj ^ 

where Nbb = (233.9 it 2.5) x 10^ is the number of BB events, i is an index that runs through 
the three Dg decay modes, Cj is the acceptance in the i-th mode, and BRj is the product of all 
secondary and tertiary branching fractions (sec Tabic 1). 

We account for systematic uncertainties by numerically convolving C{N) with a Gaussian dis- 
tribution with width determined by the total systematic uncertainties (Table 2) in the two modes, 
including the 1.1% uncertainty in Nbb added in quadrature. We find limits B{B~ Dg(f)) < 
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1.8 X 10"*^ and B{B- D*f(p) < 1.1 x 10"^ at the 90% confidence level. As in Section 3, 
these limits are calculated using B{Ds </'7r) = (4.8 ± 0.6)% from Reference [7]. If we were 
to use the value B{Ds (frir) = (3.6 it 0.9)% from the Particle Data Group, we would find 
B{B~ — > Dg(f)) < 2.6 X 10~^ and B{B~ D*g~ (f) < 1.7 x 10~^. For completeness, we also compute 
B{B- Dg<p) X B{Ds 07r) < 8.1 x 10"*^ and B{B- D*g-(t>) x B{Ds -> ^tt) < 5.3 x 10"^ also 
at the 90% confidence level. 

In summary, we have searched for B Dg (f), and we have found no evidence for these 
decays. Our limits are about two orders of magnitude lower than the previous results, but are still 
one order of magnitude higher than the Standard Model expectation. Our limits, however, are 
much lower than expectations based on i?-parity violating supersymmetric models [3]. The upper 
limit in the B~ — > Dgcf) mode is also about a factor of four below the expectation from a two Higgs 
doublet model [3]. 
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Figure 2: Distribution of uies for (a) — > D^cp and (b) ^*s 4' candidates. The superim- 

posed curves are the result of the fits described in the text. The dashed red curve is the background 
contribution and the sohd blue curve is the sum of the signal and backgound components. 
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